How Stress Can Reduce Dissipation in Glasses 
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We propose that stress can decrease the internal friction of amorphous solids, either by increasing 
the potential barriers of defects, thus reducing their tunneling and thermal activation that produce 
loss, or by decreasing the coupling between defects and phonons. This stress can be from impuri- 
ties, atomic bonding contraints, or externally applied stress. Externally applied stress also reduces 
mechanical loss through dissipation dilution. Our results are consistent with the experiments, and 
predict that stress could substantially reduce dielectric loss and increase the thermal conductivity. 

PACS numbers: 63.50.Lm,62.40.+i,65.60.+a 



I. INTRODUCTION 

At low temperatures between 0.1 K and 10 K, a wide 
variety of amorphous solids exhibit a universal plateau 
in their mechanical dissipation Q~ x ~ 10~ 4 — 10~ 3 [1, 
2]. However, there are exceptions such as in amorphous 
silicon where doping with 1 at.% of hydrogen reduces 
the low temperature internal friction plateau by about a 
factor of 200 [U . In addition the dissipation in high stress 
silicon nitride (SisN4) thin films, which show no long 
range order in X-ray diffraction and TEM images, is 2 to 

3 orders of magnitude lower than in amorphous Si02 from 

4 K up to room temperature Q. Such a large effect is 
surprising since the stress of 1.2 GPa corresponds to only 
about 70 K/ atom. Even the dissipation of stress relieved 
SiaN4 has a Q that is about an order of magnitude lower 
than typical amorphous solids (ij. 

So far no theoretical explanation for these results has 
been presented. In this paper we propose that all these 
reductions in dissipation are due to stress but cannot 
be explained by one physical effect. Impurities, dopants, 
and internal bond constraints can produce internal stress. 
Externally applied stress can reduce dissipation through 
dissipation dilution [f| as Saulson has pointed out [g. 
In addition we propose that stress, whether internal or 
external, can reduce the dissipation produced by micro- 
scopic defects known as two level systems (TLS), either 
by increasing TLS barrier heights or by decreasing the 
coupling between phonons and TLS. Our goal is to urge 
experimentalists to make further measurements to quan- 
tify the role of dissipation dilution as well as to differen- 
tiate between these two possible effects of stress on TLS. 

In dissipation dilution 0] externally applied stress in- 
creases the stiffness of materials without increasing their 
loss, resulting in a higher Q. A simple example of dis- 
sipation dilution would be the increase in Q of a mass 
suspended from a lossy spring when a stiffer lossless 
spring is added in parallel to the original spring. Since 
Q = f j (A/) where f a is the resonant frequency and A/ 
is the line width (full width half max), f a , and hence 
Q, increase without increasing the damping. In the ap- 
pendix we estimate that a thin film square resonator of 



high stress silicon nitride could have a Q up to 40,000 
times higher than a hypothetical stress-relieved silicon 
nitride resonator due to dissipation dilution. This far 
exceeds the experimental factor of order 150 by which 
external stress increases Q in high stress silicon nitride 
The full enhancement of 40,000 is not realized proba- 
bly due to external sources of dissipation, e.g., clamping 
losses. 

Dissipation dilution only plays a role when there is 
externally applied stress. So even though dissipation di- 
lution can have a dramatic effect, it cannot explain why 
dissipation is lowered by an order of magnitude or more 
in materials which have no externally applied stress, e.g., 
in silicon doped with 1 at.% hydrogen Q or in stress 
relieved S13N4 Also, in addition to dissipation di- 
lution, external stress may reduce the internal friction 
arising from microscopic defects. To understand this, 
we note that Q^ 1 = A<ft where tf> is internal friction, 
and A is due to dissipation dilution and is a function 
of macroscopic parameters, e.g., elastic moduli Q. (For 
the rest of the paper, except where noted otherwise, we 
will focus on the internal friction and set A = 1 so that 
we can use Q _1 in place of 4> in order to be consistent 
with the accepted notation in the field of glasses at low 
temperatures.) We propose two possible ways in which 
stress could reduce internal friction: either by increas- 
ing the barrier heights of microscopic fluctuating defects 
or by decreasing the coupling (deformation potential 7) 
between phonons and two level systems (TLS). We fit ex- 
isting, but incomplete, experimental data on dissipation, 
specific heat, and thermal conductivity for silicon nitride 
and SiC>2, finding somewhat better fits to the Q _1 data 
of SiaN4 at high temperatures with the barrier height 
model. Further measurements could distinguish between 
these two models. 

In glasses at low temperatures, acoustic loss at low 
frequencies is attributed to TLS [2J, I8l4l2lj . While the mi- 
croscopic nature of TLS is a mystery, one can think of a 
TLS as an atom or group of atoms in a double well po- 
tential that can sit in either well. At low temperatures, 
the lowest 2 energy levels dominate. The TLS density of 
states is assumed to be uniform at energies below a few 
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Kelvin, so if stress merely shifts the density of states, 
there should be no effect. At low frequencies and tem- 
peratures, the primary mode of attenuation is relaxation 
in which the phonon at the measurement frequency mod- 
ulates the TLS energy level spacing The measure- 
ment frequency is not related to the TLS energy because 
the incident phonon can modulate TLS with any energy 
splitting. Attenuation occurs when the TLS population 
readjusts to the equilibrium Boltzmann distribution with 
the aid of the entire thermal distribution of phonons. 

Low acoustic loss could have important implications 
for dielectric loss since the two are completely analogous 
within the TLS model [H| ■ TLS with electric dipole mo- 
ments can produce dielectric loss by attenuating photons. 
So we would expect stressed dielectrics to also have low 
dielectric loss that could make them useful substrates 
to reduce loss and noise in superconducting qubit cir- 
cuits [lsl ]. For example, hydrogenating amorphous silicon 
nitride decreases its dielectric loss tangent by approxi- 
mately a factor of 50 16j . 

At low temperatures tunneling dominates, but at 
higher temperatures thermal activation over energy bar- 
riers becomes important. One possibility is that stress 
increases the potential energy barriers V which reduces 
tunneling and thermal activation, thus effectively reduc- 
ing the number of defects and the internal friction. We 
will show that this approach is quantitatively consistent 
with measurements of Q _1 in stress relieved SisN4, and, 
even if we ignore dissipation dilution and demand that 
the entire reduction be due to a reduction in internal 
friction, with measurements of in high stress S13N4. 
We use a single set of parameters to calculate the 
specific heat C(T), and the thermal conductivity k(T) 
in Si02 and silicon nitride. Since low dissipation im- 
plies a long phonon mean free path and a high thermal 
conductivity, we predict that the thermal conductivity of 
stress relieved SisN4 is an order of magnitude higher than 
amorphous Si02 from 4 K up to room temperature, and, 
if there is no dissipation dilution, the thermal conductiv- 
ity of high stress S13N4 could be even higher, potentially 
making silicon nitride a useful substrate for integrated 
circuits where cooling is important. 

The paper is organized as follows. We describe our 
calculations of the dissipation, thermal conductivity, spe- 
cific heat, and dielectric loss in section II. In section III, 
we explain our procedure for determining the parameters 
for fitting the experimental data. The results of those 
fits to the specific heat, thermal conductivity, and dissi- 
pation are presented in section IV. We discuss why the 
dissipation of stress relieved S13N4 is lower than ordinary 
materials in section V. We discuss the possibility that 
stress could reduce dielectric loss in section VI. In Section 
VII we present an alternative model for how stress could 
lower the dissipation, namely by reducing the coupling 
between TLS and phonons. We summarize our work in 
section VIII. 



II. CALCULATIONS OF DISSIPATION, 
THERMAL CONDUCTIVITY, SPECIFIC HEAT, 
AND DIELECTRIC LOSS 



Let us briefly review the TLS model [8|, |9|. The 
TLS Hamiltonian is H = H Q + H e where H a = 
(1/2) [Aa z — A D a x ] and H e — ^ea z where A is the energy 
asymmetry between the potential energy wells, A Q is the 
tunneling matrix element, 7 is the deformation potential, 
e is the strain field, and a x and a z are Pauli matrices. 
The energy eigenvalues of H a are E — ± y/ A 2 + A 2 . We 
follow Tielburger et al. 



17| and approximate the double 



well by two overlapping harmonic oscillator wells, each 
with energy level spacing Ml . The tunneling matrix el- 
ement A Q is given by the WKB approximation [l7j |: 



A 



Ml 



^ (VaTI + VK) exp(- VA 2 + A) (1) 



where A = 2V/(Ml ) and V is the height of the en- 
ergy barrier. Fits to the low temperature thermal con- 
ductivity find the TLS density of states P that couples 
to phonons to be approximately constant. However, an 
excess of local vibrational states, referred to as the bo- 
son peak, is evident at higher temperatures and energies 
ll!, 18]. We model these modes by Einstein oscillators 
with a step function in the density of states that starts 
at an energy E Q typically between 10 and 40 K (l9| . 

According to the TLS model, at low frequencies {v < 
1 THz) and low temperatures (0.1 K < T < 10 KL Q" 1 
is a temperature independent constant given by |l0| : 



Qo 1 



irPj 2 
2pv 2 



(2) 



where p is the mass density, and v is the sound veloc- 
ity. The sources of attenuation are TLS relaxation pro- 
cesses (Q^iTLs)' resonan t scattering of phonons from 
TLS (Q^sTLs) ana - Einstein oscillators (Q^ Q ) in which 
the phonon energy matches the energy level spacing, and 
Rayleigh scattering (Q^ y ) from small scatterers of size a 

such that ka ~ 1 where k is the phonon wavevector [l9j |. 
Yu and Freeman [19j found that a = = Hv/E a is 
consistently ~ 25% larger than the size [2(| of a molecu- 
lar unit for SiC>2, GeC>2, polystyrene, and PMMA (poly- 
methylmctahcrylate) . Just as in their work, we cut off 
Rayleigh scattering at E a . We include thermal activation 
as well as direct phonon relaxation in the TLS relaxation 
processes [17], and assume that the relaxation attenua- 
tion from Einstein oscillators is negligible 19]. Thus we 
can write 17, 19] 
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The attenuation due to TLS relaxation is given by 



' rcl.TLS 



2Q0 1 
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(4) 

where f VA = J^ x dV f* V dAP(A, V)/P with Knax = 
V + 6o- - -P(A, V) is the TLS distribution of A and V. 
We assume that A has a uniform distribution and V has a 
Gaussian distribution with an average Vo and a variance 
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(V-Vo) 2 
2al 



(5) 



The TLS relaxation rate r _1 is the sum of the direct 
phonon relaxation rate t^ 1 in which the excited TLS 
decays to the ground state by emitting a phonon, and 



the rate r Arr of Arrhenius activation over the barrier: 
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(6) 
(7) 

(8) 



where the sum is over the longitudinal and transverse 
phonon modes and tq = 2/f2o- For SiC>2 to = 4 x 10 _12 s. 
For ujT m <C 1, Q~} TLS ~ Qq , where r m is the minimum 
relaxation time for a TLS with energy E at temperature 
T The Rayleigh and resonant phonon scattering 

terms are given by 

Bvoo 3 
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where S K is the step height in the density of states of 
the Einstein oscillators that is used to fit the thermal 
conductivity k, and B is a constant. 

Q^ 1 is measured at low frequencies of order 1 MHz. 
Estimating the order of magnitude of the various contri- 
butions at 1 MHz and 1 K using the values of the pa- 
rameters in Table 1 for Si02 (transverse phonon modes), 
we find Q; c l TLS ~ Q- 1 ~ 6 x 10" 4 , Q'l^s ~ 
Q- 1 tanh(fto;/2fc B r) - 1 x 1CT 8 , and Q R ] y - 2 x 1CT 15 . 
Thus TLS relaxation dominates Q _1 at low temperatures 
and low frequencies where the plateau in Q^ 1 is given by 



Q 



plat 



exp 



2al 



(12) 



This replaces Eq. ((2J , and is obtained by plugging Eq. (j4]) 
into Eq. ([3]) and noting that the dominant contribution 
to the integral in Eq. Q is for V <C Vq due to the expo- 
nential dependence of r^ r and tJ 1 on V. The factor of 



exp[— (Vo/<7o) 2 ] in P(A,V) effectively reduces the num- 
ber of active TLS. 

The relaxation time r in Eq. (j4]) for Q~J TL s i s expo- 
nentially sensitive to the barrier height V because both 
the tunneling matrix element A Q in (see Eqs. (p} and 
(O) and the thermal activation time TA rr given by Eq. © 
depend exponentially on V. We assume that stress in- 
creases the barrier heights V, thus increasing the relax- 
ation times Td and TArr, and reducing the dissipation 



^ rcl. 
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In our model stress increases the av- 
erage barrier height Vo and decreases the variance cr 2 , in 
P(A,V). 

In order to determine the values of the parameters re- 
quired to fit Q , we need to fit the thermal conductiv- 
ity k(T) and the specific heat C(T). The equations for 
C(T) and k(T) are as follows. In glasses heat is carried 
by phonons [2l|. k(T) is given by 



C D (T,u>)v£(T,uj)du> 



(13) 



where ojd is the Debye frequency, and we approximate the 
phonon specific heat by the Debye specific heat Cd(T", oj). 
The phonon mean free path £ is related to Q by 

l(T, u) = Q(T, uj)v/u = Q(T, w)A/(2tt) (14) 

where A is the phonon wavelength. 

The specific heat C(T) has contributions from the 
phonons which we approximate with the Debye specific 
heat Cd, from TLS CtlSj an d from local modes which 
we model with Einstein oscillators Ceo fl9i | : 



where 



C{T) = C D (T) + C TLS (T) + C EO (T) (15) 



N ( T \ f XD , e x 

CD = V*fe) I dxAx J^W (16) 



Ctls = k B P I x 

V,A 



(e x + l) 2 6 



Ceo = n S c k B T j dx '^ e x _ ^2 



n k 2 B T (17) 



(18) 



where x = E/k B T, x = Ml /k B T, x D = Q D /T, N/V is 
the number density of formula units, and 0(E) is a step 
function. Q B is the Debye temperature. n a is the TLS 
density of states that contributes to the specific heat, and 
Sc is the size of the step in the density of states due to 
the Einstein oscillators that contribute to C(T). 

The dielectric loss tangent tan<5 is analogous to the 
acoustic dissipation Q . At high frequencies and low 
temperatures (~ 1 K), the dominant scattering is reso- 
nant scattering of photons by TLS in which the photon 
energy matches the TLS energy splitting. If the electro- 
magnetic intensity J is much less than the critical inten- 
sity J c , we are below saturation, and can use Eq. (fiTj) 
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with tan d replacing Q T ^ S TLS , and Q Q 1 replaced by [22| We assume a P7 2 value (such that Q 
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' o, dielectric - 



(19) 



where n e is the density of TLS with electric dipole mo- 
ments, p is the electric dipole moment, to is the angular 
frequency of the incident photons, e is the permittivity 
of the vacuum, and e r is the dielectric constant. Since 
the integral in Eq. (fTTj) is dominated by V <C Vq, we can 
make the approximation 



tan5 = QA eIectric exp 



V 2 
2a 2 



tanh 



2k B T 



(20) 



for the barrier height model. For the model, described 
in Section VII, where the stress modifies the deformation 
potential 7, Vq = 0, and 

tan S = Q^ dielectHc tanh (J^j (21) 



III. PROCEDURE FOR FITTING THE 
EXPERIMENTAL DATA 

a. Si0 2 

To fit the data for SiC>2 we follow Tielburger et al. 
[ItJ and set Vq = 0. Then by fitting the low tempera- 
ture plateau of Q~ x using Eqs. (2) and (12), we obtain 
P7 2 . Fitting Q -1 over the whole range of temperature 
yields oq. The temperature of the rise in Q _1 determines 
Ml /2. Since Vq and do are known, we can determine no 
by fitting the specific heat C(T) which then gives the 
value of P. P, 7, Vq, and ao determine the low temper- 
ature thermal conductivity k(T) without any adjustable 
parameters. We set the energy E a of the onset of the 
step in the density of states by E a = 6_d/(27t X 1.27) (lij 
where ®d is the Debye temperature. Fitting C(T) at 
higher temperatures determines the step Sc in the den- 
sity of states due to local modes (Einstein oscillators). 
The fit to k(T) at high temperatures gives the step in 
the density of states S K and the Rayleigh scattering pa- 
rameter B. Note that P < n and S K < Sc because 
not all of the degrees of freedom that contribute to the 
specific heat scatter the phonons that are responsible for 
the thermal conductivity. No one has tried before to see 
if one set of parameters can be used to fit the data for 
all these quantities. 



b. Si 3 N 4 

Fitting the data for S13N4 is complicated by the fact 
that measurements of Q~ 1 (T), k(T), and C(T) have not 
been made for the same stoichiometry of silicon nitride. 



for S13N4 and use Eq. (12) to fit the Q _1 (T) data for low 
stress Si 3 N 4 0] to obtain V and a . Assuming SiNi.i 5 
has the same values of Vq and <Jq as low stress Si3N 4 , we 
can obtain no and P by fitting the C(T) data of SIN1.15. 
By fitting the k(T) data of SiNi.15, we obtain P7 2 , and 
thus 7. Assuming SiaN4 has the same 7 as SiNi.15, we 
obtain P for SisN4. If this value is reasonable compared 
to the 7 values for SiNi.15 and a-Si02, we stop. Other- 
wise, we choose another P7 2 and repeat the above pro- 
cedure until we obtain a reasonable value of P. We then 
follow the procedure given in the previous paragraph to 
fit C(T) and k(T) for SiNi.15 at higher temperatures to 
obtain the values of £>, S K , and Sc- 



IV. RESULTS: FITS TO EXPERIMENTAL DATA 



SiQ 2 



Our fits to the data for k(T), C(T), and Q^ 1 for Si0 2 
and silicon nitride are shown in Figures [IJ [2] and [3] with 
the parameters given in Table 1. No one has tried before 
to see if one set of parameters can be used to fit the data 
for all these quantities. The fits to the Si02 data show 
that this can be done. 

Fefferman et al. 23j have reported that around 10 
mK, the acoustic dissipation of Si02 is linear in temper- 
ature. This linear temperature dependence is attributed 
to interactions between TLS 23| . In Figure 2] we show 
our fits to the data from 1231] . To obtain these fits we 
followed Fefferman et al. [23j and added a linear term 
Tint = bT(A / E) 2 0] where b is a constant to the ex- 
pression for the relaxation rates in Eq. ([6]). 6 is a con- 
stant. 



b. Si 3 N 4 

Fitting the data for SiaNj is complicated by the fact 
that measurements of Q _1 (T), k(T), and C(T) have not 
been made for the same stoichiometry of silicon nitride. 
Assuming that Q^ 1 = <f>, i.e., with no dissipation dilu- 
tion, our predictions for C(T) and n(T) for high stress 
and stress-relieved SiaN4 are shown in Figs. [1] and [21 
Around 3 K, k(T) for stress relieved Si3N4 is about an 
order of magnitude higher than for Si02, and high stress 
SiaN4 could be even higher, which is consistent with low 
dissipation and a long phonon mean free path. 

From Table 1, we see that our fits to C(T) for 
SiN1.15io.05 require surprisingly large values of n , the 
TLS density of states; n = 4.5 x 10 47 /Jm 3 for 200 nm 
thick films and uq — 1.5 x 10 48 /Jm 3 for 50 nm thick films, 
which are two and three orders of magnitude larger than 
values for amorphous Si02, respectively. This accounts 
for the high specific heat below 5 K. 
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FIG. 1. (Color online) C(T)/T 3 vs. T for amorphous S1O2 
and silicon nitride. Experimental data points are shown for 
50 nm and 200 nm thick SiNi.15 [U and Si0 2 . The Si0 2 
C(T) data are from [l, 26]. The solid lines through the points 
are theoretical fits. Our predictions where stress affects V or 
7 are indicated in the legend by (V) and (7), respectively. 
C(T)/T 3 curves for high stress and stress relieved SisN4 lie 
on top of each other for the barrier height model. Similarly 
for the 7 model. 
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FIG. 2. (Color online) k(T) vs. T for amorphous Si0 2 and 
silicon nitride. Experimental data points are shown for 50 
nm and 200 nm thick SiNi.15 and Si0 2 . The Si0 2 k(T) 
data are from [27I I2H]. The solid lines through the points are 
theoretical fits. Our predictions where stress affects V or 7 
are indicated in the legend by (V) and (7), respectively. At 
low temperatures k(T) for high stress SiaN4 is the same for 
the V and 7 models. Similarly for stress relieved SiaN4. 
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FIG. 3. (Color online) Dissipation Q vs. T for stress re- 
lieved Si3 N4 measured at 3.5387 MHz (solid squares) [j, high 
stress SisN4 measured at 1.526445 MHz (solid triangles) B, 
and amorphous Si0 2 (open circles) measured at 11.4 kHz [17j . 
Solid lines are theoretical fits using the model where stress re- 
duces barrier height. Dashed lines are our theoretical predic- 
tions associated with reducing 7. Dissipation dilution factor 
A = 1 in the theoretical curves. 



Our model fits the Q 1 data very well. At low temper- 
atures (T < 0.1 K), Q- 1 ~ T 3 , and we predict that Q will 
increase by up to an order of magnitude from 400 mK to 
100 mK in both stress relieved and high stress SisN4. To 
obtain an upper bound for the change in barrier height 
due to stress, we ignore dissipation dilution. In this case 
from Table 1 we see that the mean barrier height for 
high stress Si 3 N 4 is V = 3.05 x 10 4 K - 2.6 eV which 
is about 33% higher than V = 2.3 x 10 4 K - 2 eV for 
stress relieved SiaN4. These values are comparable to the 
bond energies of SiaN4 (29|. This increase in Vo is consis- 
tent with our hypothesis that stress increases the barrier 
heights. To see that these numbers are reasonable, note 
that the difference AVb in mean barrier height Vo due to 
stress is 7500 K. The applied stress is estimated to be 
about 70 K/atom. n x 10 K/(N/V) in Table 1 implies 
that 0.06% or 1 in 1700 atoms are fluctuating defects. If 
the stress is distributed nonuniformly so that each atom 
contributes, say, 6% of its stress to the defect, then 70 
K/atom x 1700 atoms x 6% = 7100 K - AVq. 



V. DISSIPATION OF STRESS RELIEVED SI3N4 

Why is Q^ 1 in stress relieved SiaN4 an order of mag- 
nitude lower than Si02? 

One might naively expect stiffer materials to have less 
dissipation by looking at Eq. ([2]) and noticing that stiffer 
materials will have a higher speed v of sound. This is 
certainly true if we compare SiaN4 and Si02- A measure 
of the stiffness of a material is the Young's modulus E. 
Silicon nitride has E = 300 GPa and v — 11.7 km/s, 



6 



10 




T(K) 

FIG. 4. (Color online) Acoustic dissipation Q^ 1 vs. temper- 
ature for Si02 at various frequencies. The data for Si02 at 
11.4 kHz is from [TtJ while the rest of the Si02 data is from 
[23| . The fits to the Si02 data with the linear term bT are 
given by dashed lines while the fits without the linear term 
added are shown as dotted lines. 



while SiC>2 is less stiff and has E — 66 GPa and v = 5.8 
km/s. (We use longitudinal speeds of sound.) However, 
P7 2 can vary from material to material and seems to be 
larger in stiffer materials. For example, PMMA is much 
softer than SiC>2 with a Young's modulus E between 1.8 
and 3.1 GPa. P7 2 for PMMA is about an order of mag- 
nitude smaller than the value for SiC>2 but the values 
of their low temperature dissipation Q" 1 are comparable. 
(P7 2 - 0.16 x 10 7 J/m 3 for PMMA, and 1.6 x 10 7 J/m 3 
for Si0 2 .) 

As another example, consider SiC>2 and GeC>2. Gc02 
is softer than Si02; the Young's modulus E — 45 GPa for 
Ge02, and E = 66 GPa for Si02 but the two materials 
have very comparable values of the dissipation plateau at 
1 K: Q- 1 ~ 4 x 10~ 4 for Ge0 2 and Q" 1 - 5 x 10~ 4 for 
Si0 2 0. P7 2 = 1.6 x 10 7 J/m 3 for Si0 2 is double that 
of Ge0 2 which has P7 2 = 0.86 x 10 7 J/m 3 19], while 
pv 2 ~ 37 x 10 6 J/m 3 for Si02 which is about 50% larger 
than pv 2 ~ 24 x 10 6 J/m 3 for Ge02. In short, the only 
way to determine the correct value of P and 7 is to mea- 
sure thermal conductivity, specific heat, and dissipation 
for samples of silicon nitride with the same stoichiometry. 
Stiffness alone is not enough to determine the parame- 
ters entering into the expression for the dissipation, or to 
account for the reduction in dissipation of stress relieved 
Si 3 N 4 . 

So we are still left with the question of why is the 
dissipation of Si3N 4 is an order of magnitude less than 
Si0 2 . The reason is that the atomic bonds are more con- 
strained in SisN4. The competition between degrees of 
freedom and bond constraints is the reason why some 
materials are good glass- formers and others are not [3(| ■ 



Each m-fold coordinated atom provides m/2 constraints 
from fixed bond lengths, and (2m — 3) constraints from 
fixed bond angles [30j. Since Si 3 N4 has 3 and 4-fold 
coordinated atoms, there are 5| constraints per atom 
which exceeds the 3 degrees of freedom per atom. This is 
more constrained than Si02 which has 3| constraints per 
atom. This increase in the number of constraints reduces 
the number of defects (TLS) and produces unrelieved 
stress that increases the average barrier height, thus de- 
creasing C(T) and Q , as well as increasing k(T). 



TABLE I. Parameters for Si0 2 , A(200nm thick SiNi.is), B 
(50nm thick S1N1.15), and Si 3 N 4 . 



Quantities^ 


SiOa 


A 


B 


Si 3 N4l 


p[10 J kg/m J ] 


2.2 


2.68 


2.68 


3.18 


Wl[10 3 iii/s£ 


5.8(L)3.75(T) 


11.0 


11.7 


11.17 


Qd[K\ 


342 


610 


649 


446 


E [K\ 


43 


76 


81 


56 


P[10 4 7Jm d ] 


0.16 


3 


10 


~ 0.39 


Sc 


1300 


7.0 


2.0 


7.0* 


S K 


250 


2.5 


1.5 


2.5* 


B[lCT 4;i s7m] 


1.7 x 10 4 


8 


6 


8* 


7 [eV] 


2.24(L)1.73(T) 


5.6 


5.6 


5.6* 


nn [K] 


12 


150 


150 


150/ 130 


2Q„ Vtt [lO""] (L) 


0.28 


68 


114 


2.13 


n [10 4 7Jm 3 ] 


2.1 


448 


1490 


58.3 /56.4 


n x 10K/(N/V)[10- a ] 


1.31 x W~' A 


1.7 


5.6 


0.59/0.57 


Vb [xl0 4 K] 





2.3 


2.3 


2.3/3.05 


cro [xl0 a K] 


0.445 


9 


9 


9/7.5 



a Density p and sound velocity v are from references 0, l25l | . Q d 
is calculated from p and v. 

b Parameters marked with * for Si3N 4 are estimated from 

SiNi.15 and Si02 while ones marked with ~ are estimated from 
other materials. Stress relieved and high stress values for Si3N 4 
are separated by / with stress relieved given first. 

c L(T) stands for longitudinal (transverse) components. If no 
data is available, we use i>t ~ 



VI. DIELECTRIC LOSS 

As we mentioned in the introduction, the dielectric 
loss tangent tan<5 is analogous to the acoustic dissipa- 
tion So if stress reduces it should also re- 
duce tan 6. We can estimate the effect of stress on 
tan 5 using the expression in the appendix. For Si02 
with n e p 2 = 1.46 x 10" 4 [H, V = 0, and e r = 3.9, 
O o :L ! ec t „ c exp[-V 2 /(2 C ri]J = Q- X dielectr . ic ~ 5x 10- 4 . 
For Si 3 N 4 with e r = 7 [31], and assuming n e is given 
by P in Table 1, p = 1 D, and using the values for 
V and o from Table 1, Q~ X dlelectrlc ~ 7 x lO" 5 and 
Q o ",Lie Ct „ c exp[-V 2 /(2a 2 )] ~ 3 x lO" 6 for stress re- 
lieved Si 3 N 4 and 2 x 10~ 8 for high stress Si 3 N 4 . Thus 
stress relieved Si 3 N 4 has the potential to lower the di- 
electric loss by 2 orders of magnitude, and high stress 
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SiaN4 could have dielectric loss that is up to 4 orders of 
magnitude lower than SiC>2. 

VII. ALTERNATIVE MODEL: REDUCED 
COUPLING 7 BETWEEN TLS AND PHONONS 

Our proposal that stress reduces the internal fric- 
tion by increasing barrier heights can be made quanti- 
tatively consistent with the data. However, there are 
other possible explanations. One is that stress decreases 
the TLS-phonon coupling 7, and does not change the 
barrier height distribution. Figs. [I] [5] and [3] show the 
results of this approach with Vq = 0, <jq = 9000 K, 
P = 4.3x 10 43 /Jm 3 , and 7 = 0.37 (3.96) eV for high stress 
(stress relieved) SisN4. The rest of the parameters are 
given m Table 1 for Si 3 N 4 . The Q- 1 fit to stress relieved 
SisN4 is reasonably good, but poor for high stress SiaN 4 
at high temperatures, indicating that this model does not 
work as well as our hypothesis that stress increases bar- 
rier heights if dissipation dilution plays no role. However, 
it is possible that some other set of values for the parame- 
ters could improve the fit to the dissipation of high stress 
SisN4. The predicted C(T) and k(T) resulting from de- 
creasing 7 is shown in Figs. [T] and [2] for both high stress 
and stress relieved SiaN4. Reducing 7 produces a ther- 
mal conductivity that is about the same as that of the 
barrier height model up to about 4 K and then, at high 
temperatures, is greater than that of the barrier height 
model by an order of magnitude or more. The specific 
heat associated with reducing 7 is about 2 orders of mag- 
nitude lower than that of the barrier height model at low 
temperatures. If stress reduces 7, the dielectric loss will 
be the same for high stress and stress relieved SiaN4 with 
Qo,dielectric ~ 7 x 10" 5 . The dielectric loss for Si0 2 will 
be the same as in the barrier height model. 

The way to differentiate between these models and to 
determine the role of dissipation dilution is to measure 
C(T), k(T), Q~ x and tan S for high stress and stress 
relieved SisN4, and determine consistent values of the 
parameters P, 7, Vq and <jq. If dissipation dilution is the 
sole cause of the reduction of dissipation by externally 
applied stress in high stress SiaN4, the thermal conduc- 
tivity, specific heat, and dielectric loss of high stress and 
stress relieved samples of silicon nitride should be the 
same. 



VIII. SUMMARY 

We have proposed three possible explanations for the 
reduction in dissipation due to external and internal 
stress. These explanations are dissipation dilution, stress 
increases the tunneling barrier Vo, and stress decreases 
the TLS-phonon coupling 7. We have used quantitative 
fits to show that these models are plausible. The only 



way to determine the respective roles of these effects is to 
determine the parameters experimentally by measuring 
the dissipation, thermal conductivity, and specific heat 
on samples with the same stoichiometry. 

It is perhaps useful to view our work in the context 
of the history to two level systems and glasses at low 
temperatures. The original model of two level systems 
was proposed by Anderson, Halperin, and Varma, and 
independently, by W. A. Phillips. It assumed a flat dis- 
tribution of the asymmetry energy and the tunneling 
barrier height of two level systems. This has been an 
enormously useful model for fitting the low temperature 
thermal conductivity, specific heat, dissipation, etc. Tiel- 
burger, Merz, Ehrenfels, and Hunklinger used a Gaussian 
distribution of the barrier height to extend the model to 
fit the dissipation over a broader temperature range. Yu 
and Freeman represented higher energy excitations with 
Einstein modes to fit the thermal conductivity and spe- 
cific heat at higher temperatures. Our paper moves the 
model forward one more step in two ways. First, we show 
for the first time that one set of parameters can be used to 
fit dissipation, specific heat, and thermal conductivity at 
both low and high temperatures. Second, we extend the 
model to include the effect of stress on two level systems. 
We propose that stress reduces the effective number of 
two level systems, either by increasing the tunneling bar- 
rier height or by decreasing the TLS-phonon coupling. 
As a result, stress would decrease the dissipation and di- 
electric loss as well as increases the thermal conductivity 
which could have important practical applications. Ex- 
amples include substrates for integrated circuits where 
cooling is crucial, and superconducting qubits where low 
dielectric noise is important. 
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APPENDIX 

In this appendix we show the calculations involved in 
our estimate of dissipation dilution. We also show the 
sensitivity of our fits to the values of the parameters. 
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ESTIMATION OF THE CONTRIBUTION OF 
DISSIPATION DILUTION 

In dissipation dilution Q materials made stiffer by ex- 
ternally applied stress without increasing their loss have 
a higher Q. If we write a complex Young's modulus 
E = Eo(l+i<fi) where <fi is the internal friction and E$ is a 
constant, then Q^ 1 = Ac/)(ujo) where A is due to dissipa- 
tion dilution and loq is the resonant frequency 0, l32| . 

We can estimate the contribution of dissipation dilu- 
tion to the reduction in dissipation of high stress silicon 
nitride by noting that the experimental geometry is that 
of a thin film square resonator [4|. The energy of the 
resonator consists of 3 parts: the kinetic energy K, the 
energy V s from stressing the material, and the elastic en- 
ergy V e i. In this case [5J 



.4 



Vei V e l 



V s + V e l V s 



(22) 



since, as we shall show, V e i -C V s . We can make the 
approximation 



Vel (h 



stress relieved 



Vs \ fo, high 



gh stress 



(23) 



where fo is the fundamental frequency of the resonator. 
The energy of a square thin film resonator has 3 contri- 
butions 1331 



H = K + V s + V e i 

1 r L °> r L v 2 
K = - I I p [u{x, y)] dxdy 



o Jo 



V °=2 



y 



T 



Ve, 



o Jo 



du 
dx 



du 
dy 



(24) 
(25) 

dxdy (26) 



V 



EI 



o Jo 

2, 



Q 2 U \ 2 / x 



dx 2 



d u 
dy 2 



n (d 2 u\ fd 2 u\ n/1 . / 3 2 u V 

W) +2{1 ~ v) [dxd-y) 



dxd$7) 



where u(x,y) is the displacement perpendicular to the 
x-y plane, L x = L y = L are the length of the sides, p is 
the mass per unit area, and E is the Young's modulus. 
T is the tensile force per unit length given by T = Sd 
where S is the stress. I — d 3 /12(l — v) where d is the 
thickness of the plate and v is Possion's ratio which is 
0.24 for silicon nitride [13]. The equation of motion for 
u(x,y) is 



d 2 u(x,y) 
' dH 



= -EI 



T 



8 4 d 4 
dx^ + du 4 
d 2 cf_ 
dx 2 dy 2 



dx 2 dy 2 



u(x,y) 



u(x,y) 
(28) 



By assuming a solution of the form u(x, y) = 
U exp [ik x x + ik y y — itot] , we find the dispersion rela- 
tion: 



p p 



(29) 



where k 2 = k 2 + k 2 . With the boundary condition 
u(0,y) — u{L,y) — u(x, 0) = u(x,L) — 0, the funda- 
mental mode corresponds to k x = k y = 2tt/L giving the 
fundamental resonant frequency: 

Z7T 




AE 



(2-Kd\ 2 fl^ 2 



12(l-i/ 2 )p V L 



L 



(30) 



where p = p/d. 

We can estimate the fundamental resonant frequen- 
cies of a square thin film resonator made of high stress 
and stress relieved SiaN4 using values from ref. [4|: S 
= 1.2 GPa for high stress silicon nitride, d = 30 nm, 
and L = 255 ^m. The Young's modulus E = 300 
GPa [34{, and the mass density p = 3180 kg/m 3 . We 
estimate / ~ 3.4 MHz for the high stress resonator 
compared with the experimental value of 1.526445 MHz 
[4|. For the hypothetical stress- relieved resonator with 
S = 0, we estimate /q ~ 17 kHz. This gives a ratio 
of (/o,hi g h strcss/Jastrcss relieved) ~ 200. Thus from Eqs. 
([22]) and ((23]) A ~ 2.5 X 10~ 5 , i.e., Q is enhanced up 
to a factor of 40,000 by dissipation dilution. However, 
experimentally [4| , the Q of high stress silicon nitride is 
increased by a factor of order 150 by external stress. The 
full enhancement of 40,000 is not realized probably due 
to external sources of dissipation, e.g., clamping losses. 



SENSITIVITY OF FITS TO PARAMETERS 

To show the sensitivity of our fits to the parameters 
P, 7, cr, fl, and the frequency /, we show in Figures [5l 
IH1 and [7] how the dissipation would change if we varied 
these parameters by a factor of 2 from the values that we 
quoted in Table 1 in the paper. 
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measuring frequency /. The black squares are the Si02 data 
measured at 11.4 kHz from Tielburger et al. [17J. The black 
solid line is the fit using the values in Table 1 with a — 445 
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FIG. 7. (Color online) Acoustic dissipation Q _1 vs. temper- 
ature for various values of a and and Vq. The black squares 
are the experimental data points for high stress S13N4 from 
4]. The black solid line shows the fit to the data. We have 
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show the sensitivity of the fit to the values of the parameters. 
The values of the other parameters are given in Table 1. 
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